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STATISTICAL  METHODS  FOR  THE  JOINT 
OIL  ANALYSIS  PROGRAM 


1.  INTRODUCTION 

This  report  contains  the  results  of  a  continuing  effort,  sponsored 
by  the  Kelly  AFB,  to  develop  statistical  procedures  to  improve  the  analyses 
of  the  spectrometric  wearmetal  readings  in  used  oil  samples;  the  results 
of  the  analyses  are  used  as  an  aid  in  monitoring  the  wear  condition  of 
equipment  with  oil -wetted  parts.  For  FY  81,  we  were  tasked  to  address  the 
development  of  (1)  a  statistical  procedure  for  the  daily  standardization 
check  of  the  Baird  Atomic  emission  spectrometers  (2)  formulas  and  tables 
to  convert  wearmetal  readings  from  an  emission  spectrometer  to  an  equivalent 
set  of  readings  from  an  absorption  spectrometer  (3)  a  moving  regression 
methodology  for  monitoring  the  condition  of  a  given  piece  of  equipment 
based  on  the  wearmetal  buildup  "signature"  exhibited  by  the  piece  of 
equipment  and  (4)  to  assist  the  Joint  Oil  Analysis  Program  Data  Coordinating 
Group  in  improving  the  program's  data  base  and  the  periodically  issued 
management  and  technical  reports. 

The  new  methodology  for  daily  standardization  check  of  the  atomic 
emission  spectrometers  is  presented  in  Section  1.  Whereas  the  current 
procedure  requires  all  laboratories  stocking  oil  standards  with  0,  3,  10, 

30,  50,  100  and  300  ppm  concentrations  the  new  one  requires  stocking 
either  0,  50  and  100  ppm  or  0,  50,  and  300  ppm  standards  only.  In 
addition,  the  current  procedure  only  ensures  that  the  accuracy  of  a 
spectrometer  is  within  bounds;  the  one  proposed  examines  both  the  accuracy 
and  repeatability  of  a  spectrometer. 


Regression  based  formulas  and  tables  for  the  conversion  of  measurements 

from  an  emission  spectrometer  into  an  equivalent  set  for  an  absorption 

spectrometer  and  vice  versa  are  in  Section  2.  When  an  aircraft  moves 

from  its  home  base  to  a  different  base  it  may  happen  that  the  two  bases 

have  the  two  different  types  of  spectrometers  respectively.  The  conversion 

formulas  will  be  useful  in  maintaining  continuous  oil  analysis  records. 

The  formulas  were  derived  using  used  oil  sample  records  from  the  correlation 

test  results  for  the  period  Mar  80  -  Feb  81.  The  coefficient  of  determin- 
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ation  R  measuring  the  goodness-of-fit  of  the  regression  was  about  .99  for 
all  the  elements  tested  except  copper.  This  would  indicate  that  the 
formulas  should  provide  reasonably  good  conversions,  in  general.  However, 
we  have  been  informed  by  Mr.  Richard  E.  Lee,  the  Army  representative, 
J0AP-‘''''C,  that  care  should  be  exercised  in  using  the  formulas;  if  the  size 
of  the  wearmetal  particles  exceeds  the  capability  of  the  emission  spectro¬ 
meter  one  can  get  misleading  results. 

Section  3  contains  the  new  moving  regression  methodology  for  tracking 
the  wear  condition  of  a  piece  of  equipment.  The  criteria  presently  in  use 
for  evaluating  the  oil  analysis  results  consist  of  trend  ».aDles  and  decision 
guidelines  that  are  derived  from  pooled  historical  data  on  all  equipment  of 
a  given  type.  We  believe  that  each  piece  of  equipment  will  exhibit  a 
unique  "signature"  of  wearmetal  buildup  that  can  be  exploited  to  improve 
the  decision  process.  To  this  end,  we  developed  a  moving  regression 
methodology  that  uses  the  most  recent  fixed  number  of  observations  to  assess 
if  a  new  observation  should  be  considered  normal,  suspect  or  indicates 
abnormal  wear  needing  attention.  The  procedures  have  been  applied  to 
several  cases  of  unit  failure  histories  that  led  to  a  JOAP  hit  using  the 
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current  criteria.  In  all  cases,  except  one,  the  moving  regression  detected 
an  abnormality  at  or  before  the  time  of  the  JOAP  hit. 

In  addition  to  the  foregoing,  we  met  with  the  data  base  committee 
of  the  Joint  Oil  Analysis  Program  Coordinating  Group,  and  provided 
suggestions  for  the  structure  of  the  data  base  and  the  types  of  managerial 
and  technical  reports  needed  for  a  successful  oil  analysis  program. 


2.  DAILY  STANDARDIZATION  CHECK  PROCEDURE 


2.1  The  Recommended  Procedure 

The  current  daily  standardization  check  allows  the  laboratory  operator 
to  select  three  concentrations  from  the  0,  3,  10,  3C,  50,  100,  300  ppm 
standards;  it  is  possible  that  a  laboratory  then  could  choose  3,  10  and 
30  ppm  standards,  say,  for  its  daily  check,  totally  avoiding  the  two 
operator-controlled  standardization  points  at  0  and  100  ppm.  The  current 
prut-odure  then  compares  the  three  single-burns  at  the  selected  levels  with 
tabled  limits,  apparently  the  same  for  essentially  all  elements.  If  one 
or  more  of  the  single-burns  does  not  fall  within  the  specified  limits,  two 
more  burns  are  made  of  the  offending  standard  and  the  average  of  three  burns 
then  is  compared  with  the  same  limits.  If  this  average  of  three  burns,  for 
one  or  more  standards  (and  presumably  one  or  more  elements)  still  lies 
outside  the  prescribed  limits,  the  operator  must  then  perform  a  complete 
standardization  procedure. 

The  new  recommended  daily  standardization  check  requires  two  daily 
burns  at  each  of  the  0  (offset),  50  and  100  ppm  levels,  (or  0  (offset), 

50  and  300  ppm  levels),  which  then  are  used  to  check  both  the  accuracy  and 
the  precision  of  the  instrument  for  all  elements  which  the  laboratory 
expects  to  monitor  on  that  day.  For  the  two  burns  at  the  same  concentration 
(for  each  element  of  interest)  the  operator  must  then 

(a)  Average  the  two  values  and  subtract  the  known  true  concentration 
(50  for  the  0  offset  value  and  for  50  ppm,  100  for  the  100  ppm,  300  for  the 
300  ppm  level).  The  magnitude  d.|  of  this  difference  (ionoring  whether  it 
is  positive  or  negative)  then  is  compared  with  the  appropriate  entry  of 
Table  1.  This  is  done  for  each  element  the  lab  expects  to  analyze  that  day. 
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So  long  as  each  computed  magnitude  is  no  larger  than  the  entr>  in  Table  1, 
the  lab  passes  the  accuracy  check  and  goes  on  to  the  precision  check  in  (b). 
If  one  or  more  elements  have  iifference  whose  magnitude  exceeds  the 
tabled  limit,  the  lab  must  make  two  more  burns  of  the  standard  (at  the 
same  concentration  level),  again  average  these  two  new  values,  subtract 
the  true  value  and  compare  the  magnitudes  (for  the  elements  monitored) 
with  the  Table  1  limits.  If  again  one  or  more  elements  produces  a  magnitude 
in  excess  of  the  value  in  Table  1,  the  operator  must  perform  a  complete 
standardization  procedure. 

(b)  For  the  two  burns  which  passed  the  accuracy  check,  the  operator 
must  now  take  the  difference  of  the  two  readings  for  each  element,  for 
each  concentration  (0  (offset),  50  ppm,  100  ppm)  or  0  (offset),  50  ppm, 

300  ppm).  The  magnitude  d^  of  this  difference  (again  Ignoring  whether 
it  is  positive  or  negative)  then  is  compared  with  the  appropriate  entry 
in  Table  2,  for  each  element  to  be  analyzed  that  day.  If  each  computed 
magnitude  does  not  exceed  the  tabled  value,  the  daily  standardization  test 
is  passed  and  the  machine  is  ready  to  use.  If  one  or  more  computed 
magnitudes  exceeds  the  tabled  limit,  two  additional  burns  must  be  made 
at  the  same  concentration*,  these  two  burns  must  then  be  used  to  repeat 
both  the  accuracy  check  in  (a)  and  this  precision  check.  If  either  of 
these  two  checks  are  failed  for  one  or  more  elements  of  interest  at  one 
or  more  concentrations,  the  operator  must  perform  a  complete  standardization 
procedure. 

The  following  example  illustrates  this  suggested  procedure. 

Example 

Laboratory  A  expects  to  analyze  only  for  Fe,  Ag,  Cr,  Cu  and  Mg 

on  a  given  day,  with  typical  sample  readings  between  10  and  60  ppm. 


Thus  for  his  daily  standardization  check  he  will  make  two  burns  of 

standards  at  the  0,  50  and  100  ppm  level,  setting  the  machine  to  the 

offset  position  for  the  0  ppm  standard.  Only  the  readings  for  the 
six  elements  listed  above  will  be  monitored.  Assume  the  results  of 
the  6  burns  are  as  listed  in  Table  3,  which  also  gives  the  computations 
for  the  accuracy  check  (a).  Since  for  each  of  the  six  elements,  for 
all  3  concentrations,  the  value  for  d.|  is  smaller  than  the  entry 
in  Table  1,  the  accuracy  check  is  passed. 

Table  4  repeats  the  same  basic  data  (two  burns  each  at  0,  50, 

100  ppm,  readings  for  the  same  6  elements)  and  illustrates  the  compu¬ 
tation  of  the  d2  values  for  the  precision  check.  These  d2  values 
are  compared  with  the  values  in  Table  2.  Notice  that  each  of  them  is 
no  larger  than  the  Table  2  value,  except  the  value  of  1,1  for  Cr  at 
0  ppm,  which  exceeds  the  tabled  value  of  1.  Thus  this  laboratory  must 
make  two  new  burns  at  0  ppm  (offset)  and  again  go  through  both  the 

accuracy  and  precision  checks  for  0  ppm.  The  two  new  burns  are  given 

in  Table  5,  which  also  presents  the  accuracy  and  precision  computations. 
Again  the  d.j  values  are  compared  with  Table  1  while  the  d2  values 
are  compared  with  Table  2.  Each  magnitude  now  does  not  exceed  the 
tabled  value  (although  d.|  =  1  for  Cr  equals  the  tabled  value) 
so  the  machine  has  now  passed  the  daily  standardization  check  and  the 
complete  standardization  is  not  required. 

If  each  machine  is  properly  standardized  each  day,  by  either  the 
reco’.nmended  daily  check,  or  if  necessary  the  full  standatc'-ation  procedure, 
it  should  be  in  proper  operating  condition.  It  follows  that  the  laboratories 
require  prepared  oil  standards  at  only  the  0,  50  and  100  ppm  levels  (or  0, 
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50  and  300  ppm  levels)  for  their  daily  operation;  provision  of  prepared 
standards  at  other  levels  (3,  10  and  30  ppm)  are  not  required.  If 
during  the  course  of  a  day's  work  it  is  desired  to  check  the  accuracy 
of  the  instrument  (as  may  frequently  occur)  the  daily  standardization 
check  may  be  repeated  at  any  one,  two  or  all  three  of  the  concentration 
levels  employed. 

2.2  Origin  of  Tables  1  and  2 

The  numbers  presented  in  Tables  1  and  2  are  derived  from  a  mixture 
of  theoretical  model  constructions  of  how  the  AE-3  readings  should  behave, 
both  from  day  to  day  and  from  one  Instrument  to  another,  and  from  an 
examination  of  actual  AE-3  readings  from  several  instruments  burning 
the  same  samples  on  two  different  days. 

The  0  (offset),  CO  ppm  and  100  ppm  readings  were  taken  from 
data  collected  by  LCDR  D.  C.  Hatcher  for  his  thesis  "Accuracy  and 
Repeatability  Indices  for  JOAP  Data",  MS  in  Operations  Research, 

September  1979,  Naval  Postgraduate  School,  Monterey.  Ca.  Twenty  six 
different  laboratories  participated  in  LCDR  Hatcher's  study.  Each 
participating  laboratory  burned  the  same  oil  sample  15  times,  on  each  of 
two  days,  on  its  own  instrument,  monitoring  up  to  20  different  contaminant 
(listed  in  Table  1);  samples  at  the  0,3,  10,  50  and  100  ppm  levels 
were  employed  in  his  study.  The  300  ppm  values  used  came  from  a 
special  separate  study  in  which  each  of  10  different  laboratories 
burned  the  same  300  ppm  sample  10  times  on  two  successive  days, 
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monitoring  the  14  contaminants  iHentit'ied  in  Table  1.  Three  of  these 


laboratories  were  also  included  in  the 


done  by  Hat 'her, 


For  each  elsment-ppm  combination  ci  merest*  the  same  basic  type  of 
underlying  model  was  assumed,  frequently  called  a  components  of  variance 
model;  this  type  of  model  allows  one  to  estimate  the  contributions  of  the 
different  laboratories  (’instruments),  of  the  days  within  laboratories, 
and  of  the  inherent  instrument  variability  to  the  observed  variability  In  the 
numbers  recorded.  Specifically,  if  we  let  y^jj^  represent  the  k~ 
reading  made  by  instrument  i  on  day  j  (for  a  given  element-ppm  combin¬ 
ation,  say  iron  at  50  ppm),  we  can  assume  that 


■'  *  -'i  *  hj  ^  Sjk  • 


The  parameter  u  represents  the  "true"  concentration  level,  .\^  represents 
the  contribution  from  laboratory  1,  ■5^j  represents  the  effect  of  day  ,j  within 
laboratory  1,  and  represents  the  observation  ery'or,  tlie  inherent  in¬ 

strument  variability  over  repeated  burns  of  the  same  sample.  I'urtliermore , 

2 

we  can  assume  the  .\^  values  are  noniial,  1ndo|)endent,  mean  0,  variance 

tlie  values  are  noni)al ,  independnit,  mean  0,  variance  and  the 

2 

values  are  normal,  independent,  mean  0,  variance  a  ;  the  specification 

of  the  model  is  completed  by  assuming  all  A^'s,  mutually 

Independent.  It  then  follows  that  any  individual  y^i^  is  normal,  mean  ii , 

2  2  2 

variance  o'  +  .  The  observed  data  can  be  used  to  estimatM 

2  2  2 

u,  a  ,  0^,  0^.  The  simplest  procedure  to  use  (and  thoj  one  employed)  to  esti¬ 
mate  tlie  variance  parameters  is  to  perform  an  analysis  of  variance  and  use 
linear  functions  of  the  observed  mean  squares  to  construct  unbiased 
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2  2  2 

estimators  for  a  ,  o^,  a^;  this  procedure  is  described  in  many  textbooks, 
such  as  Statistics  in  Research,  by  Ostle  and  Mensing,  Iowa  State  University 
Press,  Third  Edition,  1975.  As  mentioned  in  this  text,  it  is  possible 
that  this  procedure  may  give  negative  estimates  for  certain  variances, 
which  are  clearly  nonsensical;  in  such  cases  it  is  prudent  to  replace  the 

negative  value  by  zero.  In  5  out  of  the  148  cases  in  which  a  negative 

value  could  occur  in  this  study,  the  procedure  produced  a  negative  variance 

estimate;  each  of  these  was  replaced  by  zero  for  the  subsequent  computations. 

2  2  2 

Granted  that  the  data  described  can  be  used  to  estimate  o  ,  a^,  o^, 
what  then  is  done  with  these  values  to  derive  reasonable  limits  for  the 
average  (or  difference)  of  two  successive  readings  of  the  same  sample  on 
the  same  instrument  (as  given  in  Tables  1  and  2)?  Wo  shall  first  discuss  the 
derivations  for  the  average  of  two  such  readings. 

For  simplicity,  assume  the  readings  are  done  at  laboratory  1  (so  i  ■  1) 

on  day  1  (so  j  ■  1)  and  the  readings  are  the  first  and  second  of  that 

day  (k  •  1 ,  k  ■  2) .  Thus  our  two  observed  readings  are  y^.|.|  and  y^^^ 

(for  the  given  element-ppm  combination).  The  average  of  the  two  is 


^iri  ^  ^112  ^  ,  "ill  ^  "112 

2 - '^1  "  '^1  "  2 


u  O  O 

which,  from  our  model,  is  normal,  mean  u,  variance  a"  +  .5^.  The 


methodology  described  in  Ostle  and  Mensing,  quoted  above,  leads  to  using 
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=  error  mean  square 
X2  =  days  within  labs  mean  square 
X^  =  lab  mean  square 


from  the  analysis  of  variance  for  the  given  eleinent-ppm  combination;  the 
constants  a,  b  and  c  depend  on  the  number  of  laboratories,  number  of  burns 
per  day,  and  number  of  days  per  laboratory  available  It  is  a  simple  con¬ 
sequence  of  our  model  that  y^^_  and  V  are  independent  random  variables, 
as  are  X^,  X2  and  X^*  It  also  follows  that  V  is  approximatoly  a  normal 


Cm  n 

random  variable  with  mean  0^  +  0^  and  thus 


(a) 


|y 


.rz 
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2.24  +  0? 


.075 


(b) 


,  2 


.975 


It  is  a  simp'e  matter  to  use  the  observed  data  to  estimate  0^;  call  the 

estimate  «  /  Oy  (oy  is  a  quadratic  function  of  the  three  mean  squares 
mentioned  above).  Then  from  (b)  it  should  be  approximately  true  that 


(c) 


0 


,\  - 


V  1.96o, 


.975 


The  independence  results  mentioned  earlier,  together  with  (a)  and  (c), 
lead  to 


P 


^  2.24 


/v  + 


1 .960, 


(.975)‘ 


,95 


Thus  wo  would  expect  about  95,'i  of  the  individual  readings  to  be  within 


2.24  /  V  +  l.OGOy  of  the  true  sample  content  (if  the  instrument  is 
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working  correctly).  Thus,  our  model -based  approach  leads  to  the  computation 


of  2.24  /v  +  1.96ay  ,  from  the  observed  data,  as  reasonable  9S%  bounds 

on  how  close  the  average  of  two  successive  burns,  by  the  same  instrument 
on  an  arbitrary  day,  should  be  to  the  actual  true  concentration  of  the 
sample  (for  a  given  element). 

A  second  quite  different  approach  may  also  be  employed  to  derive 
reasonable  95%  limits  for  the  difference  |y^^_  -  u|.  This  is  best  de¬ 
scribed  by  referring  to  one  actual  case,  say  iron  at  50  ppm.  For  this  case 
there  were  26  laboratories  (instruments)  involved,  each  of  which  analyzed 
the  same  sample  15  times  on  each  of  two  days  (30  x  26  «  780  analyses  m 
total).  Within  each  day  we  can  pair  together  successive  analyses  (1  and  2, 
3  and  4,  ...,  13  and  14,  drop  15)  and  average  these  pairs;  thus  we  can 
compute  7  averages  for  each  lab  for  each  day,  a  total  of  7  x  2  x  26  "  364 
independent  averages  of  pairs  (for  this  element-ppm  combination).  We  con 
then  rank  these  364  averages  from  smallest  to  largest  and  ask  how  far  below 
50  ppm,  and  how  far  above  50  ppm,  must  we  go  to  include  the  346  (.95  x  364) 
central  averages.  The  distance  we  must  go  should  be  roughly  the  same  as 

2.24  /v~+' |.96ay  derived  from  our  model  and  indeed  the  agreement  is  quite 

good.  The  numbers  given  in  Table  1  were  derived  by  looking  at  botfi  tfiese 
approaches,  for  each  element-ppm  combination. 

r.rt  entirely  analogous  procedure  was  used  for  deriving  limits  for 
tiie  difference  of  two  successive  readings,  given  in  Table  2.  From  our 
model,  using  the  notation  discussed  earlier. 
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whicn  is  nonnal  with  mean  0.  variance  2o  (both  the  lab  effect  and  the  day 
effect,  as  well  as  the  true  concentration  u,  cancel  in  forming  the  dif- 
ference).  The  estimaff'  of  2a  ,  for  any  element-ppm  combination,  is 
W  =  2X^ ,  two  times  the  error  mean  square  in  the  analysis  of  variance 
mentioned  earlier.  Again,  W  and  y.j.|.|  -  y^^2  independent  and 


111 
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P(2a^  <  W/(l  -  1.96d))  =  .975, 

where  d  =  /^  ,  v  “  error  mean  square  degrees  of  freedom. 
These  two  equations,  and  independence,  lead  to 
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so  our  iiiodel  gives  2.24  / reasonable  95%  limits  for  the  dif¬ 
ference  of  two  successive  analyses  on  the  same  instrument;  these  quantities 
were  computed  for  each  element-ppm  combination.  As  with  the  averages  of 
pairs,  differences  of  successive  analyses  were  computed,  ranked,  and  used 
to  find  how  far  above  and  below  zero  one  must  go  to  bracket  95 o  of  the 
observed  differences.  Again,  the  agreement  with  the  model-based  value  of 

/  W 

2.24  /  y  was  quite  good,  resulting  in  the  values  given  in  Table  2. 

The  values  given  in  Tables  1  and  2  should  allow  a  correctly  calibrated 
instrument  to  pass  the  daily  standardization  check  about  95?i  of  the  time. 

It  follows  that  about  5%  of  the  time  the  full  standardization  procedure  will 
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u™.c.ssar(ly  be  beeulred.  ThU  proportion  (W)  con  be  becreeseb  by  in- 
creesing  the  values  presented  in  the  two  tables.  Such  an  increase  would 
also  increase  the  probability  that  an  incorrectly  calibrated  instrument, 
which  should  employ  the  full  standardization  procedure,  will  pass  the 
daily  standardization  check. 


3.  DETERMINATION  OF  CONVERSION  FACTORS  FOR  READINGS  FROM  ATOMIC  ABSORPTION 
AND  ATOMIC  EMISSION  SPECTROMETERS 


The  joint  Oil  Analysis  Program  was  established  in  1976  for  the  purpose 
of  utilizing  spectrometric  oil  analysis  to  detect  impending  failures  of  air¬ 
craft  engines,  gear  boxes,  hydraulic  systems  etc.  Presently,  there  are 
nearly  180  spectrometric  laboratories  participating  in  the  program.  Although 
a  wide  majority  of  the  laboratories  use  atomic  emission  spectrometers 
(A/E35U-3),  a  few  of  them  still  perform  the  oil  analyses  on  atomic  absorption 
spectrometers.  Because  of  the  differences  in  the  operating  principles  of  the 
two  types  uf  spectrometers  the  analytic  results  obtained  from  the  two 
machines  tend  to  be  different  even  for  the  same  oil  sample.  In  fact, 
measurements  from  an  absorption  spectrometer  tend  to  be  lower  than  those 
for  an  emission  spectrometer. 

The  spectrometric  laboratories  depend  on  weamietal  evaluation  criteria 
to  determine  the  appropriate  action  to  be  taken  based  on  the  oil  sample 
analytic  results.  The  criteria  are  presented  in  tables  of  wearmetal  ranges 
and  trends  that  are  to  be  considered  normal,  marginal,  high  or  abnormal  for 
a  given  piece  of  equipment.  The  criterion  values  are  determined  from 
historical  records  and  statistical  distributions,  based  on  a  large  number  of 
sample  analytical  results.  The  required  statistical  data  is  obtained  by 
pooling  all  available  analytic  results  on  a  given  piece  of  equipment.  In 
order  for  this  method  to  be  successful,  it  is  important  that  all  the 
laboratories  analyzing  oil  samples  for  a  particular  type  of  equipment  obtain 
similar  measurements  for  similar  samples. 

Ensuring  uniformity  in  analytic  results  for  all  laboratories  using 
emission  spectrometers  is  a  relatively  easy  problem;  proper  calibration  of 
the  spectrometer  on  a  regular  basis  should  achieve  uniformity.  However, 
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as  stated  earlier,  emission  and  absorption  spectrometers  do  not  provide 
comparable  analytic  results  and  one  approach  to  correct  this  is  to  establish 
formulas  to  convert  analytic  results  from  absorption  spectrometers  into 
comparable  results  for  emission  spectrometers.  Such  formulas  would  also 
be  useful  in  maintaining  consistent  records  (unit  failure  histories)  on 
transient  aircraft  that  move  from  one  base  to  another  and  their  spectrometers 
are  different. 

The  Naval  Weapons  Engineering  Support  Activity  (NAVWESA),  Washington  D.C. 
in  1975  demonstrated  the  feasibility  of  using  statistical  regression  tech¬ 
niques  to  determine  the  Conversion  formulas.  The  investigators  Lynch  and 
Short  [3]  used  limited  data  from  the  interlaboratory  correlation  program  to 
demonstrate  the  techniques  and  recommended  a  more  ''xtensive  study  to  establish 
the  conversion  formulas.  Cueller,  Jr  [2]  of  the  Southwest  Research  Institute 
also  studied  this  problem  in  1975  and  recommended  the  use  of  the  following 
conversion  formulas: 

A/E35U-3  =  2AA  =  y  A/E35U-1 

Whereas  the  NAVWESA  study  indicated  the  need  for  a  separate  formula  for 
different  wearmetals,  the  Southwest  Research  Institute  results  are  to  be 
applied  uniformly  to  all  wearmetals.  One  of  the  tasks  of  our  research 
project  is  a  further  investigation  of  the  problem  for  the  purpose  of 
establishing  firm  conversion  formulas.  We  extracted  the  data  for  our 
analysis  from  the  interlaboratory  correlation  test  results  for  the  period 
Mar.  1980  -  Feb.  1981.  ^e  selected  six  elements  viz.,  iron,  silver, 
aluminum,  chromium,  cooper  and  magnesium  for  this  study.  Approximately 
140  AE  laboratories  and  25  AA  laboratories  contributed  to  the  data  base. 

Each  month  the  laboratories  analyzed  two  independent  oil  samples.  The 
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results  were  summarized  by  averaging  over  all  AE  laboratories  and  over  all 
AA  laboratories,  respectively.  Thus,  the  data  for  our  analysis  consisted 
of  24  observations  (2  per  month  for  12  months)  for  the  AE  laboratories  and 
24  observations  for  the  AA  laboratories  (Table  6)  for  each  of  the  six 
elements. 

We  used  regression  analysis  to  fit  a  straight  line  to  represent  the 

relationship  between  the  AA  and  the  AE  readings  and  the  results  were 

surprisingly  good.  With  the  exception  of  copper,  a  straight  line  provided 

2 

a  very  good  model  with  coefficient  of  determination  R  values  above  .96; 

for  copper  R  was  .75.  The  results  of  the  regression  analysis  are  presented 

in  Table  7  and  Figures  1-6  are  scatter  plots  of  the  data  and  the  goodness  of 

fit  of  linear  models  is  quite  apparent  from  these  plots,  We  are  of  the 

opinion  that  the  linear  models  described  in  Table  7  can  be  used  safely  to 

convert  an  AA  or  AE  reading  into  an  equivalent  AE  or  AA  reading  for  the 

elements  iron,  silver,  aluminum,  chromium  and  magnesium.  For  copper  also 

we  suggest  the  use  of  the  linear  model  although  it  may  not  be  totally 

satisfactory  for  converting  an  AE  reading  into  an  equivalent  AA  reading. 

To  make  it  more  convenient  for  spectrometric  laboratory  personnel  to 

use  the  results  we  have  prepared  for  each  of  the  six  elements,  ready  lockup 

tables  converting  AA  data  sets  into  equivalent  AE  data  sets  and  vice  versa. 

These  are  presented  in  Tables  8  and  9. 

Under  the  interlaboratory  correlation  progran.  tho  laboratories  also 

analyze  two  synthesized  samples  each  month.  We  exam'ned  the  relationship 

between  the  measurements  on  AA  machines  and  those  froM  the  AE  machines. 

Here  again,  we  found  that  linear  models  provide  a  practically  perfect  fit 
2 

with  an  R  value  greater  than  .98  for  each  of  the  six  elements.  The  data 
and  the  regression  results  are  presented  in  Table  10  and  11. 
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TABLE  FOR  CONVERSION  OF  AA  READINGS 
INTO  EQUIVALENT  AE  READINGS  (PPM) 
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TABLE  FOR  CONVERSION  OF  AE  READINGS 
INTO  EQUIVALENT  AA  READINGS  (PPM) 
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TABLE  OF  MEANS:  SYNTHESIZED  SAMPLES 
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SYNTHESIZED  OIL  SAMPLE  MEANS 

LEAST  SQUARES  PIT  OF  A  STRAIGHT  LINE 

Y  •  a  +  aX  (X  -  AA  MEAN;  Y  ■  AE  MEAN) 
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4.  STATISTICAL  CRITERIA  FOR  THE  EVALUATION  OF  OIL  SAMPLES 


4.1.  Description  of  the  procedure. 

The  current  wearmetal  methodology  consists  of  a  wearmetal  evaluation 
criteria  and  trend  tables  as  well  as  a  decision  making  guidance  table  for 
each  type  of  equipment  included  in  the  JOAP  program  [1].  The  evaluation 
criteria  and  the  trend  tables  have  been  derived  using  historical  data  on 
each  type  of  equipment.  When  a  new  sample  is  received  and  analyzed  the 
evaluator  uses  the  evaluation  criteria  and  trend  tables  to  make  one  of 
several  available  recommendations  such  as: 

a)  Continue  routine  sampling 

b)  Submit  another  sample  to  confirm  "abnormality" 

c)  Increase  the  sampling  frequency  or 

d)  Take  a  prescribed  maintenance  action  immedia ;.ely. 

It  is  true  that  historical  data  can  be  use^  effectively  to  construct 
the  tables  described  above.  However,  we  believe  that  each  individual  piece 
of  equipment  will  exhibit  a  unique  "signature"  of  wearmetal  buildup  that 
can  also  be  exploited  to  improve  the  decision  process.  In  an  earlier 
technical  report  [2]  a  regression  based  approach  to  monitoring  the  wear¬ 
metal  buildup  rate  for  a  given  piece  of  equipment  was  discussed.  We 
now  describe  a  modified  version  of  the  procedure  that  is  capable  of 

a)  recognizing  the  trace  or  signature  of  metallic  contamination  which 
is  normal  for  a  piece  of  equipment  in  good  operating  condition  and 

b)  is  capable  of  detecting  changes  in  the  signature  that  may  indicate 
an  abnormal  condition.  The  basic  assumption  behind  the  methodology  is 
that  an  engine  in  proper  working  order  accumulates  metallic  contaminants 

in  its  lubricating  fluid  in  some  well  behaved  pattern  and  any  major 
deviation  from  the  pattern  could  be  attributable  to  an  abnormal  condition 
that  should  be  corrected. 
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Let  K  be  the  number  of  wearmetals  to  be  monitored  by  the  procedure. 

While  the  equipment  is  in  good  working  order  the  lubricating  fluid  is  sampled 
at  prescribed  intervals  (such  as  after  each  flight  of  an  aircraft)  and 
analyzed  on  a  spectrometer.  Whenever  an  oil  sample  is  taken  the  number 
of  accumulated  hours  since  the  last  oil  change  is  recorded.  The  initial 
data  for  the  "moving  regression  methodology"  consists  of  N  observations 
taken  at  times  t.|  1  t2  £  •  •  •  £  tj^j  hours  since  the  previous  oil  change.  The 
observations  are  the  spectrometer  readings  Y..,  i  =  1,2,  ....  N; 

*  J 

j  =  1,2,  ...,  K  on  the  K  wearmetals.  We  shall  assume  that 


yj(t.)  =  cj  .  ej 


t.  +  e .  . 
1  TJ 


that  is,  each  of  the  K  elements  may  change  linearly  in  concentration  with 
flight  times;  8.  represents  the  rata  of  change  in  concentration  for  element 

J 

j  and  uj  represents  the  initial  concentration  of  element  j  in  the  fluid. 

J 

The  e.  .  are  measurement  errors,  caused  by  inhomogeneities  of  the  fluid  and 

•  J 

the  inherent  randomness  of  the  spectrometer  readings  etc.  It  is  also  assumed 
that,  for  each  element  the  are  independent,  normal,  each  with  variance 

'  J 

2 


Standard  linear  regression  methodology  can  be  used  to  estimate  the 

2  2 
parameters  ,  8j  and  a.;  these  estimates  are  denoted  by  a.,  b.  and  st 

respectively.  Rather  than  examining  the  regression  results  on  the  wearmetals 

individually  we  propose  to  find  a  single  linear  combination  Y(t)  =  ^c.Y.(t) 

which  is  most  informative,  in  a  sense,  regarding  rates  of  accumulation  of 

the  contaminants  and  regarding  changes  in  rates. 
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For  any  set  of  constant  weights,  Cp  C2,  ....  our  assumed  model 
implies 

V(t,)  =  ICjVj(t,) 

=  Xcja.  *  t.Jc.Sj  +  c^e, .  . 

That  is,  we  would  expect  Y(t)  to  also  be  a  linear  function  of  the  time  of 
observation  t.  If  a  new  oil  sample  is  received  taken  T  hours  after  oil  change, 
the  value  we  would  expect  for  Y{T),  based  on  the  initial  data  set  is 

A 

Y(T)  =  Ic.a.  +  T^c.b..  The  estimated  signal-to-noise  ratio,  frequently  employed 

J  J  J  J 

in  many  technical  areas,  is,  for  a  sample  taken  T  hours  after  oil  change, 

V(t)  , 

v^Var  Y(t)  Var  Y . (t) 

J 

choosing  the  c.'s  to  maximize  this  quantity  gives  the  "most  sensitive" 

J 

linear  combination  Y(T)  in  some  sense.  It  was  shown  in  the  NPS  Technical 
Report  [2]  cited  earlier  that  the  optimal  coefficients  are 


/var”  Y.(t] 


where  Var  Y.(T)  = 

J  J 


N+l  ,  (T  - 

T  Ktr'ti! 


To  keep  the  range  of  the  Y(T)  values  on  roughly  the  same  scale  as  the  original 
readings,  it  is  convenient  to  normalize  the  weights  c.  so  that  they  add  to 

J 

one;  this  is  done  by  letting  the  modified  coefficients 
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c . 


Granted  the  weights,  C^,  Cn.  . . . ,  Cj^  have  been  determined,  they  can  be  used 
to  determine  upper  bounds  which  should  not  be  exceeded  (with  a  specified 
probability)  for  the  composite  reading  taken  T  hours  after  oil  change.  We 
propose  the  use  of  a  student's  -  t  approximation  to  the  distribution  of 
''''(T)  =  IC.Y.(T)  to  determine  the  bounds;  the  resulting  bounds  will  be  of 

J  J 

the  form 

Y(T)  +  t^  /id  Var  Yj.(t) 

4*  U 

where  is  the  upper  y'"  percentile  of  a  student's  t-distribution  with 
N-2  degrees  of  freedom. 

For  any  given  serial  number,  known  to  be  in  good  working  order,  oil 

samples  are  analyzed  at  specified  time  intervals  (such  as  after  each  flight 

of  an  aircraft)  giving  N  sets  of  contaminant  readings.  The  regression 

2 

estimates  a..,  b.,  s,  are  computed  as  are  the  weights  C, ,  C,,  ...,  C^,.  For 

*J  J  ^ 

a  new  sample  taken  T  hours  since  oil  change,  confidence  bounds 
Y.g(T),  Y.gg(T)  and  Y.gg(T)  are  computed.  These  bounds  an  the  values  which 
would  be  exceeded  only  10  percent,  5  percent,  and  1  percent  of  the  time 
respectively  If  the  equipment  is  in  good  order  and  may  be  used  for  making 
decisions  as  to  whether  the  new  sample  is  to  be  considered  unusual  or  not. 

For  example,  a  'K'  code  may  be  assigned  if  Y(T)  is  between  the  95  and  99 
percent  bounds,  a  'T'  code  is  assigned  if  Y(T)  exceeds  the  99  percent  bound, 
increase  the  sampling  frequency  if  Y(T)  is  between  the  90  and  95  percent 
bounds  and  continue  normal  sanipl’'nc  if  Y(T)  is  below  the  90  percent  bound. 

If  a  decision  is  made  rot  to  recomiiiend  a  special  maintenance  action  and  to 
continue  s.npling  the  initial  data  base  is  updated  pending  the  arrival  of  the 
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next  sample.  The  update  is  performed  by  deleting  the  oldest  observation 
corresponding  to  the  sampling  time  t-j  and  adding  the  most  recent  obser¬ 
vation  obtained  at  time  T  to  the  data  base.  Once  a  new  sample  is  analyzed 
the  regression  computations  are  reinitiated  to  determine  if  this  sample  is 
to  be  treated  as  unusual. 

If  the  decision  is  to  taxe  a  special  maintenance  action  such  as  an 
overhaul  then  we  recommend  that  a  completely  new  data  base  consisting  of 
N  observations  be  created  before  initiating  the  regression  methodology. 

It  can  happen  that  even  when  a  given  piece  of  equipment  appears  to  be  in 
good  working  order,  the  oil  gets  changed  as  part  of  routine  maintenance. 

If  fiis  happens,  the  contaminant  readings  in  samples  taken  prior  to  the 
oil  change  and  subsequent  to  the  oil  change  tend  to  be  quite  different,  the 
latter  ones  being  much  smaller.  To  achieve  compatibility  between  the  two 
observations  it  is  proposed  that  the  post  oil  change  readings  and  the  ob¬ 
servation  times  be  added  to  the  last  observation  prior  to  oil  change  and  the 
corresponding  observation  time  respectively  before  pen  'rming  the  regression 
analysis.  This  special  handling  is  to  be  repeated  with  each  of  the  first 
N-1  observations  subsequent  to  the  oil  change.  At  this  juncture  a  new  set 
of  N  observations  all  taken  subsequent  to  the  oil  change  will  be  available 
as  initial  data  and  normal  processing  can  be  resumed.  This  special  handling 
of  the  data  after  an  oil  change  can  add  a  systematic  bias  if  the  new  oil 
already  has  nonzero  levels  of  the  contaminants.  Ideally,  it  would  be  pre¬ 
ferable  to  measure  these  new  oil  contaminant  levels  and  subtract  them  from 
the  new  observations  before  initiating  the  special  treatment. 

Details  of  the  computational  steps  involved  in  applying  the  moving 
regression  methodology  are  presented  in  Section  2  and  Section  3  contains  an 
example  illustrating  the  computations  and  the  possible  decision  actions. 
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If  this  methodology  were  to  be  implemented  operationally  it  would  be 
necessary  for  the  SOAP  laboratory  personnel  to  have  access  to  some  type  of 
micro-computer.  Then,  a  computer  program  can  be  written  and  stored  which 
can  be  invoked  whenever  the  methodology  is  to  be  applied.  In  fact,  the 
computer  calculations  can  be  simplified  considerably  by  using  vector  and 
matrix  notation  thereby  handling  all  the  wearmetals  simultaneously.  Once 
an  appropriate  computer  program  has  been  developed,  all  that  the  laboratory 
personnel  need  to  do  is  input  each  new  observation  and  invoke  the  program. 

The  methodology  was  applied  to  several  sets  of  unit  failure  history 
data  (December  1979)  for  various  aircraft  and  in  all  but  one  case  the 
moving  regression  put  out  an  "alert"  at  the  same  time  or  earlier  as  the 
JOAP  laboratory  recommended  teardown  (Tables  13-21). 

4.2.  Computational  Procedures  for  the  Moving  Regression  Methodology 

The  moving  regression  methodology  provides  a  means  to  track  the 
"signature"  exhibited  by  a  sequence  of  observations  (spectrometic  results) 
of  a  piece  of  equipment  in  the  SOAP  program.  The  observed  signature  is  then 
used  to  oetermine  if  a  new  observation  is  to  be  considered  "extreme"  and 
warranting  possibly  some  maintenance  action.  The  notation  and  the  compu¬ 
tational  formula;  are  described  below. 

Notation 

K  is  the  number  of  wearmetals  under  scrutiny.  For  example,  if  the 
wearmetals  of  interest  are  iron,  silver,  aluminum,  chromium,  copper,  and 
magnesium,  then  K  =  6. 

N  is  the  number  of  observations  to  be  used  in  the  moving  regression. 

i  is  the  subscript  used  to  index  the  observations  or  samples.  If 
N  -  10  observations  are  to  be  used  in  the  moving  regression  then  i  takes 
the  values  1 ,  2,  . . . ,  10. 


'->«*  .■'il'iitnr  ii  ~  ' ■  .'^£ks 
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j  Is  the  subscript  used  to  index  the  K  wearmetals.  If  K  =  6,  then 
j  takes  the  values  1,  2,  ....  6. 

The  initial  data  consists  of  K  sets  (one  for  each  wearmetal)  of  N 
observations  Y^,  Y^,  ....  Yj^.  Tt  avoid  cumbersome  notation  the  wearmetal 
subscript  j  has  been  omitted  temporarily.  After  some  of  the  initial 
computations  are  completed  the  wearmetal  subscript  j  will  be  introduced. 
The  initial  computations  described  below  will  be  exactly  the  same  for  each 
of  the  K  wearmetals,  using  the  appropriate  wearmetal  data  set. 

4 

The  above  observations  are  obtained  at  times  t^ ,  t^.  tj^  where 
the  times  t^  are  measured  in  "hours  since  previous  oil  change".  It  is 

assumed  that  these  times  form  a  non-decreasing  sequence, 
i  .e. ,  t,  <  to  <  • . .  <  t.  .  Note  that  these  times  t,.  will  be  the  same  for 

I  -  1 

each  of  the  K  wearmetals. 

Initial  Computations 

For  each  of  the  K  wearmetals  compute  the  following: 


l  t,l  I  t  1  I  V,:  I  I  t.' 
i=i  ^  i=l  ^  i=l  ’  i=»l  ^  i  =  l  ^ 
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(3)  Regression  coefficients  and  variance  estimates; 


Slope  parameter  b  = 
Intercept  parameter  a  =»  7  - 


Variance  estimate  S 


2  „  ^22  “  hi 

e  “N  (N“-  2T“S,, 


The  wcarmotal  subscript  j  will  be  introduced  at  this  stage.  At  the  end 
of  the  above  calculations  there  will  be  available  K  sets  of  regression 
estimates  (sj*  tij,  j  »  1,2,  ....  K. 

Main  Computations 

A  new  observation  or  a  sample  taken  T  hours  since  oil  change  is  available 
and  it  is  to  be  decided  whether  this  observation  is  to  be  considered  "extreme" 
or  "abnormal".  The  new  observation  will  consist  of  the  spectrometer  readings 
on  each  of  the  K  weamietals  ,  Yy2.  •••.  Y^|^  .  Instead  of  examining  the 
individual  wearmetal  readings  separately  the  methodology  determines  an 
"optimal  linear  combination"  over  the  wearmetals  to  be  used  as  a  decision 
function.  The  necessary  computations  are  listed  below. 

(4)  Check  if  T  >  tj^  the  time  of  the  immediately  preceding  observation. 

This  should  be  true  unless  an  oil  change  occurs  subsequent  to  the  previous 
observation.  If  T  <  t^,  proceed  to  step  (5);  otherwise  go  to  step  (6). 


(5)  Define: 


T'  =  T  +  t^  and 


. 

Proceed  to  the  next  step  (6)  and  perform  the  calculations  replacing  T 
with  T'  and  with  Y|j.  This  step  artificially  creates  continuity  in  the 
data  by  assuming  that  had  not  the  oil  change  occurred  the  new  observation 
time  would  have  been  T'  and  the  new  observation  Yj^ ,  j  =  1,2,  K. 
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Each  time  an  oH  change  occurs  this  step  is  to  be  repeated  with  the 


first  N  observations  at  the  end  of  which  a  new  set  of  initial  data  all 
obtained  after  the  oil  change  will  be  available  for  normal  processing. 
(6)  Compute: 

gj  >  aj  t  Tbj 


N  <■  1  H(T  -  il 


11 


=0  ■  Sj/I'j 


and 


V  /  J, 


for  j  ■  1,2,  ....  K. 

The  Cj  are  the  coefficients  in  the  optimal  linear  combination. 

(7)  Compute  the  observed  value  of  the  optimal  linear  combination  (OLC) 

'T-j,  ■ 

the  expected  value  o^"  the  OLC 


A 


K 


i  Tbj)  . 


the  variance  of  Y  (T) 


j  =  l  ^  ^ 


and  the  upper  confidence  bounds 
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where  t  ^q,  t^gg  and  t  gg  are  the  upper  90  ,  95  and  99  percentile  of  a 

student's  t  distribution  with  N  -  2  degrees  of  freedom;  these  values  can 
either  be  looked  up  from  appropriate  tables  or  computed  directly  using  an 
appropriate  program. 

(8)  Compare  the  observed  value  of  the  OLC,  Yy  with  the  upper  confidence 
bounds  to  make  decisions  on  the  appropriate  laboratory  recommendation.  As 
an  example,  If  Yy  Is  below  the  90  percent  bound  assign  an  A  code  to  continue 
routine  sampling,  Increase  the  sampling  frequency  for  closer  monitoring  of 
the  equipment  If  Yy  Is  between  the  90  percent  and  95  percent  bounds;  If  Yy 
Is  between  the  95  percent  and  99  percent  bounds,  assign  a  'K'  code  calling 
for  a  retagged  sample  and  If  Yy  exceeds  the  99  percent  bound  assign  a  'T' 
code  recommending  Immediate  grounding  for  appropriate  maintenance  action. 

(9)  Update  the  Initial  data  for  the  moving  regression  as  follows: 

If  the  current  observation  did  not  lead  to  a  special  maintenance  action 
add  the  current  data  to  the  Initial  data  base  deleting  the  oldest  records. 
Thus,  the  sampling  times  which  were  t.| ,  t^,  ....  t^  would  become 

tg,  ....  t^^,  T.  Each  wearmetal  concentration  Is  treated  the  same  way; 
for  any  given  wearmetal  (suppressing  the  j  subscript)  the  records  which  were 
Yp  Y^,  ...,  Y|^  become  Y^,  Y^,  ....  Y^^,  Yj^^^  where  Y|^^^  is  the  latest 
observed  spectrometer  reading  for  the  wearmetal. 

If  the  current  observation  Is  sufficiently  unusual  to  warrant  a  special 
maintenance  action,  as  stated  in  Section  1,  accumulate  a  new  data  base 
consisting  of  N  observations  all  taken  subsequent  to  the  maintenance  action 
before  Invoking  the  moving  regression. 


4.3.  An  Example 

An  example  illustrating  the  computational  steps  is  presented  below. 

The  example  starts  with  an  initial  data  set  and  each  time  a  new  observation 
is  made  available  the  moving  regression  computations  are  initiated.  The 
data  for  the  example  (Table  1)  uses  the  unit  failure  history  of  a  T-38B 
aircraft  (serial  number  400287)  for  the  month  of  December  1979.  The  results 
of  each  of  the  computational  steps  will  be  shown  for  the  first  iteration 
only.  After  that  only  the  results  of  the  most  significant  steps  will  be 
presented.  Six  wearmetals  viz.,  iron,  silver,  aluminum,  chromium,  copper, 
and  magnesium  will  be  tracked  in  the  example.  It  should  be  pointed  out 
that  in  this  example  all  the  aluminum  readings  are  Identically  zero. 

However,  we  did  Include  this  wearmetal  since  it  is  one  of  the  more  commonly 
occurring  wearmetal  in  the  SOAP  program  and  also  to  Illustrate  that  its 
inclusion  has  no  effect  on  the  analysis. 

Number  of  samples  in  the  moving  regression  N  =  5 

Number  of  wearmetals  K  =  6 


Time  since 
oil  change 


86 

106 

128 

205 

223 


INITIAL  DATA _ 

Spectrometer  Readings 


Fe 


6 

5 

7 

9 

10 


Ag 

1 

0 

1 

0 

0 


A1 

Cr 

Cu 

’i3 

^4 

^•5 

0 

1 

2 

0 

1 

2 

0 

2 

3 

0 

2 

6 

0 

3 

6 

Unit  Failure  History  for  December  1979 
Aircraft  Type:  T-37B  Serial  Number:  400287 


Time  Since 

Oil  Change 

Fe 

Ag 

A1 

Cr 

Cu 

Mg 

86 

6 

1 

0 

1 

2 

2 

106 

5 

0 

0 

1 

2 

1 

128 

7 

1 

0 

2 

3 

2 

205 

9 

0 

0 

2 

6 

2 

* 

223 

10 

0 

0 

3 

6 

3 

244 

9 

1 

0 

2 

7 

3 

264 

11 

1 

0 

4 

9 

4 

ii-k 

1 

4 

0 

0 

1 

2 

1 

24 

13 

0 

0 

1 

2 

1 

44 

23 

0 

0 

2 

3 

3 

64 

27 

0 

0 

2 

3 

2 

86 

32 

0 

0 

2 

3 

2 

105 

37 

0 

0 

2 

4 

2 

•kick 

108 

40 

0 

0 

3 

4 

2 

* 

The  first  five  rows 

are  used 

as  the 

initial 

data 

base 

in  the  example 

k* 

It  appears  that  an 

oil  change  occurred  just 

prior 

to 

this  observation 

kkk 

The  unit  failure  history  indicates  a 

JOAP  recommended 

teardown  and 

the  auxiliary  drive 

FT  inlet 

was  found  to  be  worn 

out 

, 

Table  12 
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Computational  Results 


Fe 

Cr 

Cu 

Step  1 

748 

748 

748 

748 

748 

748 

= 

126,770  126 

.770  126 

,770 

126,770  126,770  126,770 

Zv, 

= 

37 

2 

0 

9 

19 

10 

zv? 

B 

291 

2 

0 

19 

89 

22 

Zt,v, 

= 

6017 

214 

0 

1527 

3336 

1613 

Step  2 

1 

- 

149.6 

149.6 

149.6  149.6 

149.6 

149.6 

? 

= 

7.40 

.40 

0 

1.8 

3.8 

2.0 

= 

74346  74346  74346 

74346 

74346  74346 

^22 

s 

86 

6 

0 

14 

84 

10 

= 

2409 

426 

0 

903 

2468 

585 

Step  3 

2.55 

1.26 

0 

-0.02 

-1.17 

0.82 

= 

0.032 

-0.006 

0 

0.012  0.033 

0.088 

S2. 

ej 

= 

0.775 

0.398 

0 

0.357  0.077 

0.202 

A  new  observation 
the  data  is: 

taken 

at  T  =  244  hours  since 

oil 

change  is 

available 

and 

= 

9 

1 

0 

2 

7 

3 
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step  4 


Step  6 


Step  7 


Step  8 


Step  9 


Fe 

bl 

Cu 

Since  T  =  244  > 

ts  =  223 

we  proceed 

to 

step  6 

10.459 

-0.141 

0 

2.947 

6.934 

2.743 

h. 

J 

0.953 

0.427 

0 

0.364 

0.249 

0.647 

10.97 

-0.33 

0 

8.10 

27.85 

4.24 

^3  = 

0.216 

-0.006 

0 

0.159 

0.548 

0.083 

Observed  value  of  the  optional  linear  combination  (OLC)  is  6.341 
Expected  value  of  OLC  =  6.756 

The  90,  95,  and  99  percent  confidence  bounds  are  7.353,  7.614, 
8.405 


The  observed  value  of  OLC  does  not  exceed  any  of  the  three  bounds; 
no  maintenance  recommendation  is  made. 


Update  the  initial  data  base  to 


^1 


Y 


i2 


Y.„ 

i3 


106  5 

128  7 

205  9 

223  10 

244  9 


0 

1 

0 

0 

1 


0 

0 

0 

0 

0 
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Note  that  the  first  row  of  the  original  initial  data  base  has  been 
deleted  and  the  last  row  of  the  new  data  base  is  the  latest  observation. 


The  next  observation  is 


Fe 

Cr 

M2. 

T  =  264 

11 

1  0 

4 

9 

4 

a . 

J 

= 

2.481 

0.282  0 

0.562 

-1 .631 

-0.010 

b. 

J 

0.030 

0.001  0 

0.008 

0.035 

0.012 

C. 

J 

0.098 

0.008  0 

0.054 

0.724 

0.115 

Observed 

value  of  OLC 

=  8.286 

Expected  value  of  OLT. 

=  7.517 

Confidence  bounds  are 

7.639,  7.849,  8.487 

The 

observed 

OLC  exceeds 

the  90  and  95  percent  bounds 

but  not 

the  99  percent 

bounds; 

assign  a 

'K'  code  calling  for 

a  retagged  sample  (the  data 

for  this  examp 1 

has  been  taken  from  the  December  1979 

unit  failure 

history 

records 

;  hence  no 

retagged  sample  would  have  been  submitted.)  Since 

no  special  maintenance 

action 

has  been 

recommended  the  new  initial  data  would  consist  of 

Vl 

Y.„ 

i3 

!i4 

Y-r 

i5 

128 

7 

1 

0 

2 

3 

2 

205 

9 

0 

0 

2 

6 

2 

223 

10 

0 

0 

3 

6 

3 

244 

9 

1 

0 

2 

7 

3 

264 

11 

1 

0 

4 

9 

4 
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The  next  observation  and  the  associated  computations  are  as  follows: 
T  =  1  Y, .  =  4  0  0  1  2  1 

T  =  1  implies  that  an  oil  change  occurred  subsequent  to  the  previous 
observation.  Hence,  we  proceed  to  step  5  and  let 

T'  =  264  +  1  Yj.  =  n+4  1+0  0+0  4+1  9+2  4  +  1 

=265  =15  1  0  511  5 

The  application  of  the  moving  regression  methodology  results  in 

C.  =  0.349  0.025  0  0.078  0.354  0.194 

J 

Observed  OLC  =  10.521 
Expected  OLC  =  7.557 

Confidence  bounds  8.261  8.569  9.503 


The  observed  OLC  exceeds  all  three  bounds;  the  appropriate  recommeridation 
would  be  to  ground  the  equipment  for  appropriate  special  maintenance  action. 

In  the  actual  unit  failure  history  data  used  in  the  example, routine  sampling 
was  continued.  We  will  therefore  assume  that  a  judgemental  decision  not  to 
take  special  maintenance  action  has  been  made.  The  data  base  will  be  updated 
by  adding  the  data  T'  and  and  deleting  the  data  corresponding  to  the 
sampling  time  t  =  128. 

The  next  4  observations  are; 


T  =  24 

13 
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T  =  64 

27 
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3 

2 

T  =  86 

32 

0 

0 

2 

3 

2 
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Proceeding  through  step  5  the  regression  results  will  be; 


T  Observed  OLC  _ Confidence  Bounds 


24 

5.257 

5.735 

6.036 

6.950 

44 

7.251 

7.860 

8.250 

9.434 

64 

10.727 

13.825 

14.370 

16.022 

86 

12.221 

14.948 

15.439 

16.929 

In  each  one  of  these  ses  the  observed  value  of  the  OLC  does  not  exceed 
any  of  the  confidence  bounds.  At  this  stage,  five  observations  all  obtained 
subsequent  to  the  oil  change  are  available  and  (until  another  oil  change 
occurs!  step  5  should  be  bypassed  in  future  computations,  with  the  new 
initial  data  base  as: 
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The  next  two  observations  are: 


T  = 

105 

37 

0 

0 

2 

4 

2 

T  -- 

108  V,.  = 

40 

0 

0 

3 

4 

2 

and  the  regression  results  are 


T  Observed  OLC  Confidence  Bounds 


105 

6.821 

7.874 

8.203 

9.202 

108 

10.1G6 

10.460 

10.836 

11.974 
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The  unit  failure  history  indicates  that  at  thfs  point  the  SOAP  laDoratory 
recommended  a  teardown  and  it  was  observed  that  the  auxiliary  drive  FT  inlet 
was  worn  out  resulting  in  a  "JOAP  hit".  Although  the  regression  methodology 
did  not  identify  the  final  sample  as  unusual  the  method  d-id  put  out  an 
alert  earlier  at  T  -■  1. 
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SERIAL  NUMBER;  103975 


Table  19 


AIRCRAFT  TYPE:  B-52G  SERIAL  NUMBER:  627986 


/  Oil  Change 

//  JOAP  HIT  -  Aux  Dr  GRBX  Ft  Inlet 
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AIRCRAFT  TYPE;  B-526  SERIAL  NUMBER:  627986 
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